Analytic Investigation of the Branch Cut of the Green Function in Schwarzschild 

Space-time 



The retarded Green function for linear field perturbations in Schwarzschild black hole space- 
time possesses a branch cut in the complex- frequency plane. This branch cut has remained largely 
unexplored: only asymptotic analyses either for small-frequency (yielding the known tail decay at 
late times of an initial perturbation of the black hole) or for large-frequency (quasinormal modes 
close to the branch cut in this regime have been linked to quantum properties of black holes) have 
been carried out in the literature. The regime along the cut inaccessible to these asymptotic analyses 
has so far remained essentially unreachable. We present a new method for the analytic calculation of 
the branch cut directly on the cut for general-spin fields in Schwarzschild space-time. This method 
is valid for any values of the frequency on the cut and so it provides analytic access to the whole 
branch cut for the first time. We calculate the modes along the cut and investigate their properties 
and connection with quasinormal modes. We also investigate the contribution from these branch 
cut modes to the self-force acting on a point particle on a Schwarzschild background space-time. 



The study of spin-field perturbations of black holes is important for many reasons. Classically they are important, 
for example, for investigating the classical stability of black holes, for the detection of field waves emitted by black 
holes and for the calculation of the self-force on a point particle moving in a black hole background space-time (which 
serves to model a black hole inspiral in the extreme mass ratio). Black hole perturbations are also important for 
understanding the quantum properties of black holes. 

A crucial object for the study of Schwarzschild black hole perturbations is the retarded Green function of the wave 
equations they obey. These equations may be separated by performing a Fourier transform in time and a multipole 
decomposition in the angular separation of the spacetime points. Thus, the calculation of black hole perturbations is 
reduced to that of the Fourier modes in the complex-frequency (w) plane followed by a sum/integral of the modes. 
Leaver [2 [5] deformed the Fourier integral along the real- frequency axis into the complex- frequency plane, thus picking 
up the singularities of the Fourier modes of the Green function. These modes possess two types of singularities as 
functions of complex frequency: an infinite number of simple poles (the so-called quasinormal mode frequencies) and 
a branch cut (BC) which lies on the negative- imaginary axis (NIA). Leaver showed that the two main contributions 
to the Green function then come from a series of modes (quasinormal modes, QNMs) at the poles and an integral of 
modes around the branch cut, which we shall refer to as BC modes. While QNMs have been extensively studied (see, 
e.g., [3] for a review), very httle is so far known about the BC modes. 

To date, only the leading asymptotic behaviour of the BC modes for small frequencies along the NIA has been 
studied at length in the literature. This small-frequency regime in the BC is known to yield a leading power-law tail 
decay at late times of an initial black hole perturbation (see, e.g., the pioneering work by Price [HE], details of the tail 
at large radius in [T] and details at arbitrary radius as well as a higher-order logarithmic behaviour in [SlI^). The BC 
modes for large frequencies along the cut have only been studied by Maassen van den Brink and by the authors [3] ■ 
In [2j it was shown that the BC modes at large frequencies lead to a divergence in the BC contribution to both the 
Green function at 'very early' times as well as to the black hole response to a noncompact Gaussian distribution as 
initial data (it is expected these divergences in the BC contributions are cancelled by similar divergences in the QNM 
contributions). The fact that highly-damped QNMs approach the BC enabled [SHE] to apply the large- frequency 
asymptotic analyses of the BC to the calculation of highly-damped QNMs. These modes have been associated to 
quantum properties of black holes (e.g., see [TOj in relation to black hole area quantization and llj in relation to 
Hawking radiation in the case of rotating black holes). 

To the best of our knowledge, the only investigations of the BC modes for frequencies which are neither asymptot- 
ically large nor small (we will refer to this regime as the 'mid'-frequency regime) are the following ones, which were 
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carried out in the gravitational case only. The BC in the Green function modes is due to a corresponding BC of a 
particular solution, gg{r,uj), of the radial equation, Eq.(2.4) below. In [12j[T3] the authors obtained the radial solution 
gi{r,u}) for frequencies near, but off, the NIA via a numerical integration of the radial equation. They thus calculated 
the radial solution on both sides of - but away from - the NIA, evaluated the difference and then extrapolated it onto 
the NIA, thus obtaining the BC 'strength'. This is a rather tricky numerical evaluation, since the difference in values 
of gi between the two sides of the NIA becomes exponentially-small as the frequency approaches the NIA. The only 
other investigation of BC modes in the 'mid'-frequency regime on the NIA was carried out by Maassen van den Brink 
who, in a different and impressive work (14], performed an asymptotic analysis of the BC modes about the so-called 
algebraically-special frequency ujas 



The algebraically special frequency lying within the 'mid'-frequency regime on the NIA occurs only for the case 
of field perturbations of spin s — 2 (axial gravitational) and correspondingly the rest of this paragraph applies to 
this case only. The BC modes have a distinct 'dipole-like' behaviour near ujas, unlike at other frequencies [T2l I13j. 
The algebraically special frequency, though not a QNM itself for axial gravitational perturbations (it is a QNM for 
polar gravitational perturbations) [13] is intimately linked to QNMs: the dipole-like behaviour of the BC modes may 
be explained in terms of poles in the 'unphysical' complex-frequency Riemann sheet [T21[I3|. Furthermore, a QNM 
frequency very close to (or exactly equal to) ujas marks the start of the highly-damped region of QNMs (e.g., [3]). 
As the rotation of the black hole is increased from zero (i.e., the Schwarzschild case studied in this paper), multiplets 
of QNMs emerge from - exactly at or very near to, depending on the azimuthal angular number - the algebraically 
special frequency ujas |121 1131 117| . at least in the case of the lowest multipole angular momentum number £ — s ~ 2. 



No analytic method exists so far for calculating the BC of the Green function in the 'mid'-frequency regime (except, 
as mentioned above, near ujas for s — 2). However, the above works (see Ref.18 in [T2] and Sec. VI [2]) suggest the 
tantalizing possibility of calculating the BC by expressing the BC 'strength' via a convergent series of irregular 
confluent hypergeometric functions evaluated directly on the NIA, so that no extrapolation onto the NIA would be 
required. In this paper we take up this suggestion. Thus, we provide a new method for calculating analytically the 
BC modes for general integral spin directly on the NIA for arbitrary values of the frequency. We prove that our new 
scries for the BC modes is convergent for any values of the frequency along the NIA, thus providing analytic access 
for the first time to the whole 'mid'-frequency regime. We note that our method is also valid in the small- and large- 
frequency regimes, but it is not useful there since convergence becomes slower as the frequency becomes small while, 
for large-frequencies, the BC modes grow and oscillate for fixed radii. Asymptotic analyses are therefore necessary in 
practise in these regimes. 



We calculate the BC modes using our new method and we investigate their properties and connection with QNMs. 
We also re-analyze the so-called Jaffe series (which is a series representation of the radial solution which is purely 
ingoing into the event horizon and possesses no BC) and, in particular, the behaviour of the Jaffe coefficients. Finally, 
we apply our calculation of the BC modes to investigate their contribution to the self- force (see, e.g., [IB]) acting on a 
point particle moving on a Schwarzschild background space-time. In [B] we 'sketched out' the main idea for our new 
method for the calculation of the BC modes for arbitrary frequency, in this paper we 'flesh out' the details. We note 
that the method we present here provided the results for the plots of quantities in the 'mid'-frequency regime in [9], 
where it was shown that these 'mid'-frequency results overlap with the large-frequency asymptotics presented there. 
In [7] we will present a thorough small-frequency analysis of the BC modes and we will show that these 'mid'-frequency 
results also overlap with that analysis in the small-frequency regime. 



In Sec|TT]we introduce the main perturbation equations and expressions for the Green function modes. In Se c|III| 
we present the various series representations which we use for the calculation of the BC modes; in particular, Eq.(3.6) 
is the new series that we derive and use for the calculation of the pivotal quantity, the BC 'strength'. In Sec IV 



we 



analyse the so-called Jaffe coefficients a„ (in particular, we correct the large-n asymptotics of these coefficients given 



in the literature) , which are fundamental in the calculation of all the series representations we use. In Sees V ■ VIII 



we calculate the various quantities required for the BC modes and these modes themselves. In Sec IX we investigate 
the contribution of the BC modes to the self-force. Finally, in Appendix [A| we give some properties of the irregular 
confluent hypergeometric function, which we require for the calculation of the BC modes in the main body of the 
paper. 



In this paper we take units c = G = 1. We will frequently use a bar over a quantity to indicate that it has been 
made dimensionless via the introduction of an appropriate factor of the radius of the event horizon, r/j = 2M, where 
M is the mass of the black hole. 
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II. BRANCH CUT 



After a Fourier-mode decomposition in time t and a multipole-^ decomposition in the angular distance 7, the 
retarded Green funcction for hnear field perturbations in Schwarzschild space-time is expressed as 

1 ^ 1 /"OO+ZC 

Gret{x,x') = ^^(2£+l)P,(cos7)G'r*(r,r';i), GT\r,r'-t) = — doj Ge{r,r';u;)e-'^\ (2.1) 

^ ' ^ J -oo+ic 



where c > and the Fourier modes of the Green function are given by 



Gf(r,r';c.) = 



W{uj) 



where r> = max(r, r'), r< = min(r, r') and r is the Schwarzschild radial coordinate. The function 

W{uj) = W[gi{r,uj)Ji{r,uj)] = gtf( - fi>g[, 



(2.2) 



(2.3) 



where a prime indicates a derivative with respect to r*, is the Wronskian of two solutions /^(r, w) and gi{r,Lu) of the 
following second order radial ODE: 



dr 



2 +cj2-y(r) 



m . (1 - ^) 



Ui{r,uj) = 



(2.4) 



where = 



rfi In 



1 J is the radial 'tortoise coordinate', rh = 2M is the radius of the event horizon and M is 

obey the 'physical' boundary conditions 

-> —00, (2.5) 
f* +00 (2.6) 



the mass of the Schwarzschild black hole. For w £ K, the solutions fi and 
of, respectively, purely-ingoing waves into the black hole 



and purely-outgoing waves out to radial infinity. 



ge{r,uj) e 



— ^ 00 



(2.7) 



The complex-valued coefhcients and A'^'^Jj are, respectively, incidence and reflection coefficients and it is straight 
forward to check that W — —2iujAY^^. The boundary conditions (2.51 and (2.71 also define, respectively, the radial 
solutions fi and gi unambiguously for Im(w) > when G M. In Im(w) < 0, with r.^ G M , the solution gi must be 

defined by analytic continuation. 

The parameter s = 0,1,2 in the potential in Eq.(2.4) is the helicity of the field perturbation, to which, with an 
abuse of language, we will refer to as 'spin': s = 2 corresponds to axial - also called 'odd' - gravitational perturbations 

1 to electromagnetic perturbations [20. and 

or ' 



2 corresponds to axial - 

(in which case Eq.(2.4) becomes the Regge- Wheeler equation [19]), s 

gravitational perturbations obey the Zerilli equation [3TJ [52] 



s = to scalar perturbations [HIS]. Polar ~ or "even' 
and solutions to this equation can be obtained from the solutions, and their radial derivatives, to the Regge- Wheeler 
equation ^23j. At the algebraically special frequency uias [SHIS]) however, this relationship between solutions to the 
Zerilli equation and solutions to the Regge- Wheeler equation becomes singular. In Figjljwe plot the potential V{r) 
for some token values of spin s and multipole number £. 

It can be shown 24J that the radial solution fi has no branch cut in the complex-frequency plane whereas ge has a 
branch cut down the negative imaginary axis NIA (see, e.g., Eq.( |3.2[ ) below). This BC in g£ can be explained [25l [26] 
in terms of the radial potential (minus the centrifugal barrier) falling off slower than exponentially at radial infinity. 
On the other hand, the exponential decay with of the potential near the horizon leads to a series of poles in fi down 
in the NIA (see Secjv] below). We note that the Wronskian W 'inherits' the BC from gi and the poles in the NIA 
from fg. We define AA{h') = A+{—ii') — A-{—iiy) for any function A — A{uj) possessing a BC along the NIA, where 
A±{-'ii') = lim£^o+ A{±e — iiy), with > Q. We will equally refer to both quantities lu and v = iu a.s 'frequencies'; 
we note that > along the NIA. 

We note the symmetries 



giir,Lo)=gUr,-io*), Mr,uj) ^ {r, ~lu*), W{lj) ^ W* {-oj*) if G 



(2.8) 
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FIG. 1. Radial potential Eq.(2.4| as a function of r,/M. Continuous blue curve: 
1=1; dashed red curve: s 
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2, £ = 2; dot-dashed green curve: 



which follow from the radial ODE (2.4 1 and the boundary conditions (2.5) and (2.7 1. These symmetries lead to 
gi- = and W+ = Wt if e M, so that the branch cuts of gt{r,ijj) and W along the NIA are only in their 
imaginary parts, their real parts having no branch cut. In particular, then, the absolute value of the Wronskian, \W\^ 
has no BC. 



The Fourier integral along the real frequency line in Eq.( 2.1 ) can be deformed into the complex- frequency plane PQH]. 
The two main contributions to G™* are, then, a series over the residues at the poles of the Fourier modes G^(r, r'; uS) 
(the QNM frequencies, which are located at the zeros of the Wronskian W) and an integral around the BC. The 
branch cut contribution G^'~^ to the retarded Green function is given by 



G 



BC 



oo 

(r, r'; tA) = — E(2^ + l)^£(cos7)Gf ^(r, r'; t), Gf ^(r, /; t) 



where the BC modes AGf can be expressed as [SI [T^ 



AGi{r,r';iy) = -2wft{r,-w)fi>[r',~ii') 



2^ 



dv AGe{r,r';iy)e~ 



We denote the function q{i') as the branch cut 'strength' as it is defined via the equation 

Agi{r, v) = iq{iy)ge{r, +w). 



(2.9) 



(2.10) 



(2.11) 



where, here, Agg{r,v) = gg^{r, —iv) — gt^{r,—iv) (the extra tilde in the notation is justified in the next section). 



From the symmetries (2.8) and the fact that gi{r, (that is, gi evaluated on the positive-imaginary axis) is real- 



valued it follows that q[i7j is also a real-valued quantity. We note that '— ig' here corresponds to the quantity 'iiT' in 
Eq.31 1 . 

As mentioned in the Introduction, we will use a bar over a quantity to indicate that it has been made dimensionless 
via the introduction of an appropriate factor r^, e.g., f = r/r/j, w = wr/j, D = vr^, etc. 



III. SERIES REPRESENTATIONS FOR THE RADIAL SOLUTIONS 



If one tried to find the radial solution fi or gi in the region Im(a;) < by naively solving numerically the radial 



Eq.(2.4) and imposing the 'boundary conditions' (2.5) and (2.7), respectively, one would run into computational prob- 



lems. The reason is that these 'boundary conditions' are exponentially dominant over the other, linearly independent 
solution at the radial endpoint where the condition is imposed (that is, at — >■ —oo for ff and at — >■ -l-oo for gi). 
Therefore, if one tried to numerically integrate the radial ODE starting with the 'boundary condition' at one endpoint 
towards the other endpoint, any accidental inclusion - no matter how small - of the other, wrong solution would grow 
exponentially and so would the numerical error. There are various methods around this problem. For example, one 
could solve the radial equation in the region Im(a;) > 0, where the boundary conditions are well-posed, and then 
analytically continue onto the region lm{ijj) < 0. Also, Leaver's Eqs. 32-36 [T provides a framework for calculating 
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the BC contribution to the retarded Green function. However, this method is rather difficult to implement (except 
in the asymptotic small-i/ regime) due to the presence of Leaver's 'phase parameter', which is required because of the 
use of a particular series representation for gj^ in terms of Coulomb wave functions. In this paper we choose to use 
certain series representations for and gi which do not involve Leaver's 'phase parameter' and which we show are 
convergent in the desired region on the frequency plane. 

Leaver j2l] provides various series representations for the radial solutions and g^. All calculations of the BC 
modes in this paper are carried using a specific choice of series representation for each one of the two solutions. 



which we give in Sees III A and IIIB However, while the new series representation for Agi (and therefore for the BC 



'strength' q{v)) which we present in Sec IIIC is fundamentally based on our choice of series representation for g£, our 



calculation of the BC modes does not depend in an important way on the specific choice of series for calculating ff. 
one could just as well use any different method valid in the 'mid'-frequency regime for calculating fi. We present the 
various series that we use in the following subsections and we investigate their convergence properties in the following 
sections. 



A. Series for fe 



In order to calculate the radial function fi, we will use the well-known Jaffe series 

Mr,u:) = (r- l)-^'"f2-e'-V,(r,w), 

oo / 1 " 

Ji{r,uj) = ^ J,i/(w), Jn,e = a„(cj) I 1 - - 



(3.1) 



We will refer to the complex- valued coefficients a,i(aj) as the Jaffe series coefficients, even though they also appear 
as coefficients in the series representation that we will use for g^, Eq.(3.2) below. The Jaffe series coefficients are 



functions of the series index n, the frequency uj and, although not indicated explicitly, the multipole number £ and 
the spin value s. The Jaffe series coefficients satisfy a 3-term recurrence relation which we give and analyze in the 
following section. The initial value oq remains undetermined by the recurrence relation; the specific value oq — e^^*'^ 
yields the desired normalization (2.5) for fi and, therefore, this will always be our choice of value for ao when using 
the Jaffe series for fg. 



B. Series for ge 



Our choice of series representation for ge is also given in [24] : 



9i 



(r,i:j) = r^+* (r - 1) e^'^''he{r,u!), hi{r,u!) = ^hn,i{^), hn,e = an{uj)Tj^ 



n=0 



T^-^ = {-2iLU + l)„ Uis + l-2iuj + n, 2s + 1, -2iujr), 



where a„(cj) satisfy the same recurrence relations as the Jaffe series coefficients a„ in Eq.(3.1 1 but it is ao 7^ ao 
is, a„ and a„ only differ by an overall normalization factor which we give below. The series (3.2) has been broadly 



(3.2) 



that 



ignored in the literature, possibly due to the fact that the irregular confluent hypergeometric [/-functions are rather 
hard to manage. We will refer to Eq.(3.2) as the 'Leaver-f/ series'. 

It is clear from the Leaver-[/ series Eq. p.2[ ) and the properties of the irregular confluent hypergeometric function [57] 
that the radial solution gi{r,uj) has a branch cut running along the line cjr : — > —00 ■ i. If r > 0, then gi{r,uj) has a 
branch cut along the NIA, w : — > —00 • i. 

The principal branch of U{a,b,z) is given by arg(z) e ( — tt, + 7r]. Therefore, we can evaluate directly on the NLA 

the confluent hypergeometric [/-function appearing in Eq.(3.2) and calculate the corresponding T^"-* via Eq.(3.2|. 



That is, ri ^ may be evaluated on the NIA and its value will correspond to the principal branch value, i.e., to the 
limiting value as the frequency oj approaches the NIA from the third quadrant in the complex-frequency plane. The 
corresponding value of gg will then give g£-^. provided that the series Eq.(3.2) converges. It will be understood, when 



we do not say it explicitly, that any quantities possessing a BC along the NIA which are evaluated on the NIA via 
the use of Eq.(3.2| will correspond to their limiting value approaching the NIA from the third quadrant. 

In order to check what boundary condition the Leaver-U series (3.2) satisfies for r 00, we use Eq.13.5.2 [281 
we obtain 



and 



(r,w) 



{f - ly"^ {-2iLur)- 



-l+2ia~ A 
floe 



for 



00 and \tt/2 — arg(a;) — arg(r)| < 37r/2 (3.3) 
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Therefore, when r > 0, the Leaver-t/ series yields the asymptotics 

gi,{r,uj) {~2iLuy~^+^''^hoe+''^''' for r ^ oo and |7r/2 - arg(a;)| < 37r/2, r>0 



(3.4) 



We note that Eq.(3.4) does not agree with Eq.75 [21]; we beheve that Eq.75 :25 is missing the first factor on the right 
hand side of Eq.(|3.4[). 

With the specific normalization choice of clq = (— 2ia;)+*+^~^*" the function calculated using the Leaver-[/ 



series satisfies the desired normalization Eq.(2.7l; therefore, this will always be our choice (different from the choice 
oo = e~^"^ above for the Jaffe series for fi) when using the Leaver-[/ series for gi. We note that a„ themselves have 
a branch cut along the NIA (this was already noted in Ref.18 of [12]), as we have 



Aa„ = - 1] a„_ = [e^^^ + l] a„+ 

where a„± = lim(;_j.o+ a„(a; = ±e — iv) and where we have assumed s G Z. 



(3.5) 



C. Series for Agi 



From Eqs.(l3.2|) and (|A4|) it follows that 

OO 

Age{r, v) = ge+{r, -iv) - ge^{r, -iv) = f^+' (f - l)""" e'^''A/i^(r, v), A/i^ (r, v) ^K, 



n=0 



_ 2^ie-2^'-e"(-+i-2g) ^ (-l)"r(l + n-2i?)C/(s-n + 2P,2s + l,2H 



r(l - 2P) 



r(l + s + n - 2P)r(l ~ s + n-2i>) 



(3.6) 



This is a series for calculating A^^ by evaluating q uant ities directly on the NIA. The principal branch is to be taken 
for the confluent hypergeometric {/-function in Eq.(3.6|. 



D. Series for Agi 



The series in this subsection, which we denote by /S.gi, would correspond to /S.gi if the coefficients a„ did not have a 
branch cut; specifically, we may view Ag^ as the discontinuity of gi across the NIA if we replace a„ by a„ in Eq.(3.2 1. 



Since that is not actually the case, we will not be using the series for l^gi anywhere. However, we inc lude it here for 
completeness, as the factors Ti"*"' in the terms of this series satisfy the same recurrence relation (Eq.(6.1) below) as 



the fact ors Tn ' and Tn^ introduced above for gi and Ag^ respectively. The solution T^'^ to t he r ecurrence relation 



Eq. (]6T|) is linearly independent from the solutions Tn ■* and r^^-* . If we replace a„ by a„ in Eq. ([3.2| , we can calculate 



the discontinuity across the NIA of the resulting quantity as: 



'■A/i^(r,i^), A/i^(r,z^) = ^ A/i„ 



^9i{r, v) = ge+{r,~ii')"^ - gi,^{r,-w)°^ = r^+'' (r - 1)" 

_ (-l)2«2^m„ (+) (+) _ r(l-2P + n)Af(l-2i> + n + s,2s + l,-2j/r) 



(3.7) 



r(i + 2s)r(i - 2v) 



r(l -2i? + n~ s) 



where we have used Eq.13.1.6 in [28]. We note the appearance of the regular confluent hypergeom etric (Kummer) 
function M in (3.7) for Agi, as opposed to the irregular confluent hypergeometric function U in Eq.(3.6) for Agi. 



E. Series for the radial derivatives 



An expression for calculating the r.^ -derivative of the radial solution fi follows straightforwardly from Eq.(3.1) 

-2 2) 



dJi 
dr 



[r — 1) r e 



— -I- 
dr f(f — 1) 



(3.8) 



— ^(n + l)a„+i 



n=0 
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In order to obtain an expression for the -derivative of gi we use Eqs.4.22-4.24 



dr^, 

dhe 
dr 



= 1- 



r ^ 

n=0 



1 + s V \ Qidhi 

+ 5f + 7- 

r r — r/i ) hi ar 



(3.9) 



Or, 



dr 



djV_ 
dr 



1-29 + n 



-2D + n+l-s , ) _ 
-2D + n + l "+i 



IV. JAFFE SERIES COEFFICIENTS 



appearing in the Leaver- L/ series Eq.(3.2) for gg satisfy the folfowing 3-term recurrence relation, 



Both the series coefficients a„(w) appearing in the Jaffe series Eq.(3.I) for fi and the series coefficients a„(a;) 

(4.1) 



a„a„+i -f /3„a„ -f 7„a„_i = 0, n = I,2, 



with a„ = for n < and where 



(n + l){n - 2D + 1) 

,2 



an 

Pn = -[2n^ + (2 - SD)n + SD^ - AD + t{l + I) + I - s~ 
7„ = — ADn + 4j>^ — 



(4.2) 



We note that although in this section we use the notation a„ to indicate a solution of Eq.(4.I) the results in this 



section apply equally to the coefficients a„ since these results are independent of the specific choice of the n = 
coefficient. 



A. Singularities of a„ 



From Eq.(4.I) it follows that, in principle, the coefficients a„ will have a simple pole where a„_i = 0, i.e., at 
n — 2v = Q. Therefore, \i v — k/2 for some fc S N then a„ will have a simple pole Vn > k (see, e.g., App.B [T2]). 
However, such a pole will not occur if at the same time it happens that Pk-ia^-i + jk-iak-2 = 0. This occurs for 
s = 2 at the algebraically-special frequency Cjas = —ii^AS, where Das = (^ ~ 1)^(^ + 1)(^ + 2)/6 [TJj. Therefore, the 
coefficients a„ do not have a pole at D = Das ior s ~ 2 while they do have a simple pole there for s = 0, 1. 

Suppose that and {c„} are two sets of solutions to a recurrence relation, then, if lim„_j.oo ^n/cn = it is 

said that 6„ are minimal and c„ are dominant. If the solution one seeks is dominant, then one can find the desired 
solution by solving the recurrence relation using standard forward recursion. However, if one wants to obtain a 
minimal solution, using forward recursion would be unstable and one must resort to finding the desired solution using, 
e.g.. Miller's algorithm of backward recursion (see, e.g., |30]). In order to investigate whether the solutions to the 



recurrence relation Eq.(4.I ) are minimal, dominant or neither, we require the large-n behaviour of the coefficients a„. 
We also require the large-n behaviour of a„ in order to study the convergence properties of any series involving these 
coefficients. 



B. Large-n asymptotics 



In order to obtain the large-rt asymptotics of the coefficients a„ we follow App.B [3T]. We thus express the asymptotic 
behaviour as the so-called Birkhojf series 



e 



Monlnn+Mi«^CogEj=oO + i"''"'^''+0("''-<'+'''''') (43) 
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for a certain chosen value of / G N, where /3 G [0, 1), Ci 7^ 0, p €E Z, p > 1 and /ij, € C for all j. Substituting this 



expression into the recurrence relation (4.1) we obtain 



48^27^ 



[315 - 1280£'' - 2560^^ + 160^^ ^ 1440^ + (7680^^ + 7680£ + 13728) D~ 



C3 

C4 = 

i 

30720^2^3/2 
(-30720^^ - 30720£ - 3200) - 10240/?^] 

2^^jQ-2 [-81 + 768£^ + 1536^3 + 480^^ - 288^ + (-1536^2 _ 1536^ - 288) D + {61Uf + 6144^ + 1152) + 
(24576s^ - 16384) + 2A576D^ + 16384^^] 

The coefficients Q are real for j even and they are purely imaginary for j odd . We note that the spin dependence 
does not appear until the term Cg. The coefficient C2 corresponds to an undetermined overall normalization and the 
' lb' s ign corresponds to the two linearly independent solutions of the recurrence relation. Since the recurrence relation 



(4.1 ) is unchanged under n — >■ ne , one solution can be obtained from the other under this change; this is essentially 



equivalent to changing the sign of Cj+i for j even in (4.4 1. On the NIA, where D > 0, the two solutions behave 
similarly (that is, no solution is dominant over the other) and an appropriate linear combination of them should be 
taken. Off the NIA, if w is not a QNM frequency then the a„ are dominant and they are generated by forward 
recursion; whereas if a; is a QNM frequency then the a„ are minimal and they can be generated by Miller's algorithm 
of backward recursion. Indeed, requiring for the solutions a„ to be minimal has become a widely used, successful 
method for finding QNM frequencies of black holes [32] . 



Finally, we note that the leading order of Eq.(4.4 1 differs from Eq.46 [24] in having a power of n equal to ^—ly — 3/4' 



instead of '— 2i/ — 3/4'; we have checked numerically for specific values of the parameters (both for w on and off the 



NIA) that Eq.(4.4) gives the correct asymptotic behaviour. 



C. Plots 



In Figj2]we show that the large-n asymptotics given in Eq.(4.4) match the exact solution to the recurrence relation 
(4.1). We note the appearance of a 'pulse', after which the values of a„ decay rapidly. Fig 3] is a 3D-plot of a„ as a 
function of both n and Mv. We have only included plots for s = 2 as representative of the behaviour of the coefficients 
a„, as the behaviour is similar for other spins. The behaviour is also similar at the algebraically-special frequency 
UJAS- At the poles described in Sec IV A the behaviour of 'sin(27ri?)a„' is also similar, except that the first n < 2D 
terms are exactly zero. 



V. CALCULATION OF fe{r,uj) 



We calculate the radial solution fe using the Jaffe series Eq.(3.1 1. As shown by Leaver in Sec. IV. A [24], the Jaffe se- 
ries is absolutely convergent Vcj e C and for any r g [r^, 00), since then: Hm„_j.oo |a„+i(l — l/f)"+-^/ (a„(l — l/f)")| = 
|1 — l/f| < 1. By the same argument the Jaffe series is uniformly convergent on r G [rh,rmax] for any finite Vmax but 
will generally not be so at radial infinity, provided the coefficients are not singular (see below). However, as shown 
by Leaver, the Jaffe series is uniformly convergent - including radial infinity - if J2n is finite, which is guaranteed 
if the sequence {a„} is minimal and this occurs at the QNM frequencies. At these frequencies it is = and 

Aout ^ — 2iu) 



As shown in Sec. IV A a„ have simple poles \fn > k when v — k/2 for some fc e N. The exception is the case 
= I' AS = i^AS for s = 2, which is not a pole. These poles carry over to fi so that this radial solution has simple 
poles dX V = k/2 (these poles of fi were shown in [331 [33] using a different method, namely, a Born series), except at 
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FIG. 2. Jalfe coefficient a„ as a function of n for s = 2, ^ = 2. Figs. (a) and (b) are for v = 0.2 and figs.(c) and (d) for v — 15.4. 



Figs, (b) and (d) are log-plot versions of (a) and (c) respectively. Blue dots: exact solution a„ to the recurrence relation (4.1 1 
with ao — 1; the continuous blue curve is an interpolation of the blue dots. Dashed green curve: large-n asymptotics Eq.(4.4l 



where we have taken a linear combination of the two linearly independent asymptotic solutions such that the linear combination 
matches the exact value of a„ at both n — 1000 and n = 2000 for v = 0.2, and at n = 145 and n = 150 for i? = 15.4. We note 
the 'pulse' centered around n = 25 in the case i? — 15.4 (the equivalent 'pulse' in the case i? — 0.2 is centered around its first 
peak at n « 20); the coefficient a„ reaches its maximum magnitude at the 'pulse' and then the magnitude decays rapidly with 



VAS when s — 2. However, the BC modes AGi are independent of the normalization of fe, and so it is useful to define 



ao = — ag sin(27riaj) 



(5.1) 



with ao = e~^"^. We denote the corresponding quantities a„, Jn,e, Ji, fe and W obtained using this normalization 
by dn, Jn/, Je, ji and W respectively. We note that at the pole v — the first nonzero value of a„ will be for 
n = k. Therefore, at v — fc/2 it is fi ^ e+''^''* as — > — oo and so fi{r,iS) (x fe{r, — w) 14j. In the particular case of 
the algebraically special frequency PaSj exact solutions to the radial equation have been found |16| . 

Therefore, as a function of w e C, the radial solution fe only has singularities at the simple poles (D — —ik/2, k £N, 
on the NIA (except at 9as for s — 2) while fg is analytic in the whole frequency plane. 

In Figj4] we illustrate, for frequencies on the positive-imaginary axis (PIA) of the complex-frequency plane, the 
convergence properties of the Jaffe series and we plot (there is no need to calculate fe on the PIA since fe has no 
poles there) as a function of \D\. In Figjsjwe do similarly but for fe on the NIA instead of fe on the PIA. In this case 
we do not plot the partial term J„ since the behaviour is essentially the same as that of a„ in Figj2|d). The radial 
derivative of fe as a function of the frequency has a similar behaviour to that of fg. In Figj6]we plot, on the NIA, 
In \ fg\ and In \dfe/dr^\ as functions of the radius: for some values of P the solution fg has a zero and for other values 
of P it does not. We note that in Figs|4j|6]we only include plots for s = 2 as the behaviour for other spins is very 
similar. In we show that the Jaffe series for fe agrees well with a small-i> series expansion. 
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FIG. 4. Construction of fi{r,bj) on the PIA using Eq.(3.1l ior s = 2, £ = 2, r = lOM. Plots (a)-(c) are with ao = 1, 



for the value lu = 9.8002z and with n on the x-axis. (a ) Partial term J„/. (b) Partial sum X^n'^o "^n.'/- ("-) Partial error 
\Jn,e/ YlZ'=o Jri',i\- (d) Log-plot of fe{r,uj) using Eq.(3.1l (blue dots interpolated by the continuous blue curve) as a function 



of and log-plot of (e'""^* +e (dashed green curve), which corresponds to the asymptotics of Eq.(2.6l ignoring the 

incidence and reflection coefficients. 
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FIG. 5. Construction of fc on the NIA using Eq. js. l| for s = 2, £ — 2. (a) Partial error ] Jn.f/ X]"'=o "^"',^1 ^ function 



of n with ao = 1 for u — 9.8002 and r — lOM. (b) from Eq.(3.1l as a function of f. The continuous blue and dotted red 
curves are calculated using the Jaffe series Eq.(3.1l for r* = 0.2M and r, — QAM respectively. Note that they are both zero 
aX u = vas because of the fact that, for s = 2, vas is not a pole of fi. The dashed green and dot-dashed brown curves are the 
large-!/ asymptotics in [9] for r« = Q.2M and r, = QAM respectively. 
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FIG. 6. Plots of In 1/^1 as_a function of r*/A'/ for s = 2, £ = 2 and (a) v = 4.4002, (b) v = 9.8002. Continuous blue curve: 

^ g-2c^r, j Dotted red curve: ln|ci/f/dr»| obtained using Eq.( 3.8 1. 
The curves for In I e'"''* + e~ 



ln|/f| obtained using Eq.( 3.1 1. Dashed green curve: ln|e'' 



Dot-dashed brown curve: In 



and its radial derivative only serve as very 



crude approximations to the large-r asymptotics of /£ in Eq. |2.5l. 



VI. CALCULATION OF m AND THE BC 'STRENGTH' 



A. Recurrence relation 



The terms in the series Eq.(3.2| for 5^, Eq.(3.6l for A(7£ and Eq.(3.7l for Ag^ all consist on the coefficient a„ times 

, Tn^ and respectively (times a factor independent of n). All three factors, 



a factor, which we denote by T, 

satisfy the same, following recurrence relation p7 1 [2^ [55] : 



(n + 1 - 2v)(n - 2v)Tn-i - (n + 1 - 2v){2n + 1 - 4P - 2vr)Tn + (n + 1 - 2i^ + s){n + 1 - 2D - s)T„+i = (6.1) 



-(-) 



andT^. 



where T„ denotes any of r„ 

In the subsections below we will show the following. On the NIA, all three factors, 

-1/4 



(-) 



the same leading-order behaviour 0{n 



^) for large-n - see, respectively, E qs.( 6.2 1, (6.31 and (6.5| below 



The solution T,[ ■* and either the solution T^^-* or T, [^^ o f the recurrence relation Eq.(|6.1 



are linearly independent . 

We show that, when solving the recurrence relation (6.1 1 for w anywhere except on the NIA, Tn ^ is a subdominant 



solution and T^^-* and TX^' are dominant solutions. In this case of w ^ NIA, if one wishes to find Tn ' , one expects that 
forward recurrence will be unstable in that the dominant solution will 'creep in' as n is increased. One should then 
use instead Miller's algorithm of backward recursion On the other hand, when solving the recurrence relation 
(6.1) on the NIA there arc no dominant nor subdominant solutions, all solutions asymptoting like 0(n~^/^). In this 
case of w e NIA, there is no danger of a dominant solution 'creeping in' and so there is no need for using Miller's 
algorithm of backward recursion for finding any of the three solutions. 



.{+) 



rp{0) 
^ 71 



and Ti+\ 



have 



(-) 



B. Calculation of i 



We need a method for calculating the radial solution gi on the PIA, as required by g(i/) in Eq.(2.11), as well as on 
the NIA, as required by the Wronskian Eq.(2.3). We will calculate directly on the NIA, as well as on the PIA, 
using the Leaver-f7 series Eq.(3.2). As mentioned in Sec. IIIB the Leaver-t/ series allows us to calculate g^ directly 

i.e., gi^. The advantage of evaluating gi on the NIA 
Secondly, on the NIA there are 



on the NIA, specifically as the limit from the third quadrant 
is twofold. First, no extrapolating procedure e — )■ O"*" onto the NIA is then needed 
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no dominant /subdominant solutions to the recurrence relation (6.1) and so there is no need for Miller's algorithm of 
backward recursion. 

The confluent hypergeometric [/-functions, however, are notoriously difficult to evalu ate. We have three options 
in order to calculate the factors in the Leaver-C/ series: (1) from their definition (3.2) and using the in-built 

[/-function in Mathematica, (2) from their definition (3.2) and using the integral representation Eq.(Al) for the U- 
function, and (3) from the above recurrence relation Eq.(6.1). We note that method (1) is highly unstable, whereas 



it is a lot more stable to calculate using method (2) (see Eq.4.16 



App.B 



1. Large-n asymptotics 



We obtain the large-n behaviour of the terms in the Leaver-[/ series Eq.(3.2) in order to investigage its convergence 
properties. From Eq.(3.2) and (A3), we have 



T{l~2D){2\iy\ry]J ^/2\U\^'^ ^ ' " ^ « ^ X - aigl fc I tt, -t-7rj, r > u 

(6.2) 



agreeing with Eq.4.19 [29]. It follows from Eq.(6.2 ) that the series for hn^i is absolutely convergent everywhere except, 
maybe, on the NIA. 

The convergence properties of the series hg = on the NIA (which would yield is the same as that for 

A/i^ = J2n ^^n,e in Sec VI C 2 and as that for Ah( — ^^n,e in Sec 
section. 



VIDl 



and so we refer the reader to this latter 



2. Results 



In FigjTjwe plot as a function of the frequency on the NIA. Since — )• e'^^ G M for large r and the asymptotic 
series for gg for large r contains only real coefficients (e.g., Sec.B.l [55^), one would expect that |Im((;£_|_)| <^ \Re{gi+)\, 
especially as u increases - this is indeed what happens in Figj?] Fig[7] shows that Im(g^+) becomes round-off error at 
some stage . In Fig. 2 [F we show similar plots for 5 = 0. 

When calculating the terms in the series Eq.(|3. 2|) fo r gi+{r, ^iv) in practise using Mathematica, the first two terms 



and r;^^i on the NIA we calculate using Eq. ( |A1[ ) and a particular splitting of some expression for the [/-function 
instead of obtaining it using Mathematica's in-built HypergeometricU function. 

In Fig|8]we plot both the radial function gg+ and its radial derivative as functions of the radius. The radial derivative 
as a function of the frequency has a very similar behaviour to that of gi^ in FigjTj 



On the PIA, as noted in SecVIBl 



solving the recurrence relation Eq.(6.1) to obtain ge is unstable since it 
corresponds to the subdominant solution and the dominant solution would be 'creeping in'. One option is to obtain 
gi{r,+ii') using Eqs.(3.2) and (A3). In Figjo] we show that using the recurrence relation is unstable whereas the 
latter option does well. The following is the method we follow: when calculating the terms in the series Eq.(3.2) for 
gi{r,+ii>) in practise in Mathematica, we use a numerical evaluation of the integral representation Eq.(Al) in order 

Vn, on the PIA instead of obtaining it by solving the recurrence relation that it satisfies. 



to calculate T, 



(-) 



C. Calculation of Age 



From Eq.13.4.15 | 
Eq.(|6ll). In Fig. [lO| we plot 



and the property T{z + 1) = zr{z) we easily see that T„ 



(0) 



Agi[r, v) as a function of v obtained with Eq.(3.6). 



satisfy the recurrence relation 
Note that, when carried out in 



practise in Mathematica, it is better to numerically evaluate the integration representation Eq.(Al) instead of using 
Mathematica's in-built Hypergeometric U function in order to calculate the two initial values T^^q and Tj^2i- 



1. Zeros and singularities of Ah^^e 



Any possible zeros and singularities of the terms Ahn/ in Eq.(3.6) may only come from the a„, the [/-function 
and the four F-functions. The [/-function has no singularities other than its branch point, and the F-function has no 
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FIG. 7. Log-plots of the radial solution ge+ as a function of for s = 2, £ — 2. Top plots (a) and (b) are for r = 2.8M and 
bottom plots (c) and ( d) a re for r = 5M. Left plots (a) and (c) are e~'^^*Ke(gi +) an d right plots (b) and (d) are lm{ge+), 
all obtained using Eq.(3.2l and finding T,[~^ by solving the recurrence relation (6.11 with Tg"' and Tj"'' calculated using 
Mathematica's in-built HypergeometricU function. Cf. Fig. 2 [9^. A similar behaviour is exhibited by the radial derivative 
dgi+ /dr-t . 
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FIG. 8. Log-plots of a nd d gi+/dr-j, as functions of r^/M for s = 2, I = 2, v = 4.4 002. These radial functions are 



obtained using Eqs. 13.21 and (3.9l and finding T^"' by solving the recurrence relation (6.11 with Tg~' and Tl~' calculated 
using Mathematica's in-built HypergeometricU function, (a) Continuous blue curve; Ke(ge); dashed green curve: large-r 
asymptotics ge ~ e""^* ; dotted red curve: |Re(d(7f+/dr*)|; dot-dashed brown curve: large-r asymptotics dgi+/dr, ~ ue"^* . (b) 
Similar as (a) but the continuous blue curve and the dotted red curve here correspond to the imaginary - instead of real - part 
of ge+ and dgi+/dr, respectively. We note that for v = 9.8002 the calculation of the imaginary part using the method as in 
these plots becomes unstable for r, 2; 8M. 



(-) 



-.(-) 



zeros. We do not know analytically the possible zeros of a„ nor of the ?7-function, so we can only determine some of 
the zeros of A/i„_£, and we cannot be sure that any of the poles we might find are not actually cancelled out by zeros 
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m 



FIG. 9. Radial solution ge on the PIA for s = 2, ^ = 2, r = 5M. (a) Continuous blue curve: Tt^ in Eq.(|3.2|) as a function of n 
for uj — 4.4002j obtained using Mathematica's HypergeometricU function (a similar curve is obtained by solving the recurrence 
relation Eq.(6.1l); this curve shows how the dominant solution 'creeps in' in t he s olution of the recurrence relation. Dashed 
red curve: similar to the continuous blue curve but T,!"^ is obtained using Eq.(All. (b) gi{r,+iu) as a function of on the 



PIA obtained by numerically evaluation the numerical representation Eq.( Al I: convergence is achieved up to large values of |a;| 
whereas using the recurrence relation Eq.( |6.1[ l the series coincides with the curve plotted up to jo; | ~ 14 but ceases to converge 
after that value. 




FIG. 10. Plot of ' 



*"* A ff^ ' a s a function of u for s = 2, I = 2, r = lOM. We have obtained l\gi from Eq.(3.6l by solving 



the recurrence relation Eq.ljBJj with Tq'"' and T]"' calculated using Mathematica's in-built HypergeometricU function. We 



.{0) 



have removed from Ay^ the large-r behaviour of Eq.(2.11l: gi{r,+iu) ~ e 
approximately 1/2, agreeing with the large-;? asymptotics of q{i') in Eq.38 |9]. 



The zeros of the curve occur at ;^-steps of 



of a„ and/or the ?7- function. 



From Sec IV A we know that a„ have simple poles Vn > k if v ^ k/2 for some fc e N (with the exception of — vas 



when s — 2). The function T{z) has simple poles at z G Z . Let us distinguish two cases: 

• Case 1 - 2P ^ U 

Neither a„ nor the F-functions have any pole. So A/i„^£ has no zeros (other than any coming from a„ or the 
[/-function) and it has no poles. 

• 1 - 2j> = -fc e u 

In the numerator, a„ has simple poles (except \iv = i/^s when s = 2) at n — k+1, k+2, . . . and F(l+n— 2P) at 
n — 0,1,2, . . . , k. In the denominator, F(l + s + n — 2P)F(1 — s + n — 2D) has double poles at n = 0, 1, 2, . . . , /c — |s| 
and simple poles at n = fc — |s| -I- 1, A; — |s| -I- 2, . . . , fc -I- |s|. Also in the denominator, F(l — 2D) has a simple pole 
Vn e N. Therefore, Ahnj has no poles. Regarding the zeros (other than any coming from a„ or the [/-function). 
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iin k, Ahn,e has double zeros atn = 0, l,2....,fc— |s| and simple zeros atn = fc— |s| + l,fc— |s|+2, ...,fc+|s|; 
the term n = k is not a zero if s 7^ and it is just a simple zero if s = 0. 

In the particular case v = vas for s — 2, a„ does not have a pole for any n € N. In this case, a similar analysis 
to the one in the above paragraph indicates that A/i„_£ has a zero there. . 



2. Large-n asymptotics 



Using Eq.(A2) and Eq.5.11.3 ||27j we obtain 



T(«) ~ "'^''^^''^^ '^' (-l)"n-V^cos (^^{ 



2n — + 1) + n [ 21^ — n — - j j , n — > 00, if s > and i^r > 

(6.3) 



Finally, together with the large-n asymptotics (4.4) for a„ we can obtain the large-n asymptotics of the terms in 
the series for Ahp: 



AK j - C(i/, r)(-l)"n"''~^e±2^^' cos ( ^/4i/r (2n - 4P + 1) + tt ( 2?^ - n - ^ ) ) , n ^ 00, if s > and i^r > 



(6.4) 



(7(i/,r) 



20Fie-'"'e^^(''+i-2p)(_2P)''+i-2i>(2i/r)-^-i/^ 
r(l - 2D) 



for fixed v. The modulus of the large-n asymptotics of Ahnj in Eq.(6.4| is basically the same as that for Ahnj in 
Eq.(6.6) below. Therefore the conclusions below Eq.(6.6) regarding tl 



le convergence of the series Ahi = J2n 



apply equally to the series Ahi = J2n ^hn,, 



D. Calculation of A(;f 



From Eq. 13.4 .1 [55^ and the property T{z + 1) = zT{z) we easily find that the T^'^ satisfy the same recurrence 



relation Eq.(6.1 1 as the Tn ■* in the series for and the Tn' in the series for Ag^ - see 



(0) 



To investigate the 



convergence properties of this series, we require the large-n asymptotics of its terms. 



Large-n asymptotics 



From Eqs.13.1.27 and 13.5.14 [H] we obtain 



TT ''-e ■'''r(2s + 1) (2i/r) - ' n '''cos 



-1/4 



( n ^ 7r(4s + l)\ 

f v/(l - + 2n) 4;/r - ^ j 



n — )■ cx). 



(6.5) 



which is valid for (n — 2i/ — s) — > +C)o, s bounded and vr £ R. This agrees with Eq.4.20 [29] (except for a typo they 
have in the sign of s inside the F-function ). [Note: Eq.(6.5 1 differs from Eq.l2 in Sec. 6. 13. 2 of [37! in having an extra 
factor 1/2 and also a power '—1/4 — s' instead of a '+1/4 — s\ We have, however, checked with Mathematica for 
certain values of the parameters that Eq.( |6.5| is the correct expression.] 

From Eq.(6.5l and the large-n asymptotics (4.4) for a„ we can now obtain the large-n asymptotics of the terms in 
the series for Agf. 



Ahn I ^ C(j/, r) n - 2P + 



1 



e^^^cos r^'^+^^ -2. 



n-2v + 



2vr 



n - 2i/ > 1, 



(6.6) 



r(i - 2v) 



17 



The ratio test yields \ Ahn+i,i/ ^hnj\ — >■ 1 as n — >■ cx), and so it is inconclusive. However, we may apply the integral 
test as follows. For v > 0, the function k{n) = \C{i',r)\n'~'^~^ = \k{n)\ is positive and monotone decreasing with n 
and it satisfies 



dn\k{n)\ ^ \C{i^,r)\ 



dn 



|C(^,r)| 



1 



lim n 

n— f oo 



< OO, 



if D>0. 



(6.7) 



Therefore \k{n)\ satisfies the integral test and so the series J^n^i''^) absolutely convergent for D > 0. Since 
l^^n.f I < 1^(^)1 for sufhciently large n, from the comparison test we have that Aft,„ ^ is also absolutely convergent 
for D > 0. Indeed, in our calculations, the series has converged for arbitrary values of i>. However, while the series is 
fast convergent for large D the speed of convergence becomes slower for smaller values P. 



E. Calculation of q{i') 



We calculate the BC 'strength' g(t^ ) usi ng Eq.( 2.11 1, where we calculate Agi using the method described in Sec VI C 
and gi{r^-\-iv) on the PIA using Eq.( |3.2[ ) with v —v everywhere. We note that an alternative method, which we 
have not explored, for calculating gi(r, -\-iv) might be to use Eq.74 [51]. In Fig. 11 we plot q{v) as a function of i>. 
This figure is to be compared with Fig. 2 [12] (also Fig. 2. [I3]). 




obtained using 



FIG. 11. BC 'strength' q{v) of Eq.( |2.1l] ) as a function of v for s = 2, I — 2. Continuous blue curve: 
Eq.(2.11 1 with the value r — 5M and gi{r, +iv) obtained using Eq.( Al I. Dashed green curve: large-i> asymptotics of Eq.38 [9] 
(see Fig. 5 [9] for a better agreement for larger values of v). Dot-dashed red curve: small-;? asymptotics of Eqs.2.2 and 2.5 [12j 
(we note that the small-t' asymptotics in [12] do not work well for other spins; see [7] for better agreement for small-i? for any 
spin). Cf. Fig.2 d2 (also Fig. 2. [E]). We note that when gi{r,+iv) is obtained using the recurrence relation Eq.(6.1l instead 
of Eq.| |Al[ ), the result for q{v) is incorrect from u ~ 11 - see Fig[9|ja). 



VII. CALCULATION OF THE WRONSKIAN 



We calculate, on the NIA, the Wronskian W+ of the radial solutions fi 



- sm(2TT iu})f£ and gi^ using the methods 
described in the previous sections: the Jaffe series Eq.(3.1| for fi, Eq.(3.8) for dfe/dr^,, the Leaver-[/ series Eq.(3.2) 
for gi^ and Eq.(3.9) for dgi^/dr^,. 

In Figsjsjjs] we plot, on the NIA, the radial solutions /^(r, —iv) and ge+{r, —iv) and their radial derivatives, which 
are required for the Wronskian. 

Let us define Wi = gt+f'i ^nd W2 = fiQe+i so that W+ — W\ — W2- Figs 12 'a) and (b) show that the magnitudes 
of the two contributions Wi and W2 are very close for all except near = 0. Therefore, the computation of 
W — Wi — W2 would require the knowledge of these two contributions to very high accuracy away from this 'window' 
near = 0. We note that for the imaginary part, for > the two contributions Wi and W2 actually add up and so 
there is no computational difficulty there either. Figs 12 'c) shows that, indeed, there is a 'window' near r* = where 
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the calculation of the absolute value of the Wronskian is reliable. We note that in this 'window' it is Ini(Ty) ^ Re(M^), 
so the imaginary part dominates but, for accuracy, the real part cannot be neglected. We have checked that a similar 
'window' nea r = occurs for different values of the spin, the multipole number ^ and the frequency on the NIA. 



In Figs 



13 



13 



c), together with Figs. 7-9 [5] where 



we plot Wi, Wi and W at r = 2.8M as functions of D. Fig 
these 'mid'-frequency results are compared to large-i? asymptotics, show that the calculation of the Wronskian at this 
value of the radius yields a reliable result. In [7] we show that the 'mid'-frequency results for the Wronskian agree 
well with small-i? asymptotics. Fig 13 c) also shows - for the particular value £ = 2 - that for s = 2 the Wronskian 
W[gi+, fi] has a zero of order one at v = vas- This is as expected because of the definition 
vas is not a pole of fi for s = 2 and it agrees with [T^] . 



and the fact that v 



■ shi{2T:iuj) fi 




M|Im(W,)| 
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FIG. 12. Wronskian W+ = W[ge+, fe;ui] = Wi ~ W2 as a. function of r./M for s = 2, £ = 2 and = 4.4002. Figs. (a) and (b): 
log-plots of the absolute values of, respectively, the real and imaginary parts of M ■ W\ = Mgt+fi (continuous blue curve) and 
M ■W2= Mftg(+ (dashed red curve). Fig.(c): Log-plot of M\W\. 



VIII. CALCULATION OF BC MODES 



We obtain the branch cut modes AG£(r, r'; v) by calculating the different quantities in Eq.(2.10) using the methods 
described in the previous sections. In particular, for the calculation of the Wronskian we have evaluated the radial 
functions fi and gi at r = 2.8M while for the branch cut 'strength' q{i') we have evaluated the radial functions at 
r = 5M. In FiglM] we plot AGi as a function of D for different spins. In the spin-2 case the plot is to be compared 



with Fig.3 [H] (also Fig.3 [13]). From Figsji] [n 
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c) respectively, the radial solution fi, the BC 'strength' q and 
the absolute value of the Wronskian \W\ all have a simple zero at v = vas in the case v — vas and s — £ ^ 2. 

c) reflects. In Fig 15 we plot 
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From Eq.(2.10) it then follows that AGi has a simple zero at that frequency, as Fig 
again AGi{r, r'; v) but in this case for larger values of the radius r': the magnitude of AGi increases rapidly with the 
radius, as expected from Figjsj^b). 

The spin-2 case is quite distinct due to the algebraically special frequency Das {— 4 when £ = 2): while the branch 
cut mode AGi is zero at = VASy AGi is particularly large for frequencies near vas- This behaviour is explained, 
in the case i = 2 and — —00 and — 00, as arising from nearby 'unconventional damped modes', that is a pair 
of poles in the unphysical Riemann sheet. The imaginary part of the QNM frequencies is negative and increases in 
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FIG. 13. Wronskian W = W[ge,ft;u}] = Wi - W2 as a. function of i? for s = 2, £ = 2, r = 2.8M (r, ^ 0.97M). (a) and (b); 
log-plots of, respectively, the real and imaginary parts of both MWi (continuous blue curve) and MW2 (dashed red curve), (c) 
Plot of the real part (continuous blue curve), imaginary part (dashed red curve) and absolute value (dot-dashed green curve) 
of A/W+/(i/ — vas)- We note the zero of order one of \W+\ at the algebraically special frequency £> = i/as- 



magnitude as the overtone number n increases, so that n is an index for the speed of damping of the mode with 
time. For spin-2, QNM frequencies approach the algebraically special frequency as n is increased from the lowest 
damped mode, n = 0, until a certain value of n, say um, whose QNM frequency is very close to ljasj for ^ = 2 it is 
nM = 9. As n is further increased from um, the real part of the spin-2 QNM frequencies in the 3rd quadrant increases 
monotonically. Therefore, in a certain sense, the algebraically special frequency marks the start of the highly-damped 
asymptotic regime for QNMs. 



In Fig 16 we plot the radius-independent quantity ^\2i'q/{MW )\\ which is the branch cut mode AGi of Eq.(2.10) 



but without the fi factors. The zeros of this radius-independent quantity correspond to the zeros of Fig 16 



shows that, for the cases with s = 1 and 2, these zeros occur with a period in D close to that of the increment in the 
imaginary part of the QNM frequencies ujqnm at consecutive overtone numbers. For s — 2, (minus) the imaginary 
part of the QNM frequencies lie close to the zeros of q{i^)- For s = 1, for which the QNM frequencies approach 
the NIA particularly fast [S], (minus) the imaginary part of the QNM frequencies lie close to the maxima points of 
'\2i'q/{MW'^)\^ which are directly related to nearby zeros of W, i.e., the QNM frequencies by definition. For s — 0, 
on the other hand, the periods of the zeros of q{iy) and of Im(wQjvM) differ slightly for mid-D while, for large-P, 
'— Im(ajQArM)' tend to lie somewhere in-between the zeros of q{v) and the maxima points of '\2i'q/{MW'^)\\ For all 
spins in the large-P asymptotic regime, the separation of the zeros of g(i/) approaches 1/2, which coincides with the 
separation in the imaginary part of highly-damped QNM frequencies for consecutive overtone numbers (see, e.g., [5] 
for the asymptotic expressions). 
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FIG. 14. Branch cut mode AGe/M of Eq.( |2.10[ ) as a function of v for r, = 0.2M and = QAM. The Wronskian has been 
calculated at r = 2.8M while qiy) has been calculated with functions at r = 5Af . (a) s = 0, £ = 1. (b) s = 1, ^ = 1. (c) s = 2, 
£ = 2 - cf. Fig.3 US] (also Fig.3 [S]). 



IX. SELF-FORCE 



The motion of a (non-test) point particle moving on a background space-time deviates from geodesic motion of 
that space-time due to a self- force (see, e.g., [IB] for a review). The self- force may be calculated via an integration 
of the covariant derivative of the retarded Green function integrated over the whole past worldline of the particle. In 
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FIG. 15. Same as Fig, 



14 



i.e., AGi/M as a function of v for s = 2, I = 2 and r* = 0.2M but here it is with r'^ = r,{r = 5M) 



in the continuous blue curve and ri = r,{r = lOM) in the dashed red curve. 



particular, for a scalar charge q moving on Schwarzschild background space-time, the /.t-component of the self-force is 
given by 



Mr)^q dr' ^^Gret{z{T),z{T')) 



(9.1) 



where z{t) is the worldline of the particle and r is its proper time. In the rest of this section we will deal with 
the case of a scalar charge (s = 0) only, although the self-force in the cases of an electromagnetic charge [s = 1) 
and of a point mass (s = 2) also involve the integration of the Green function in a similar way. We will investigate 
the contribution to the scalar self-force from a single BC multipole mode Gf-^ {r,r' \t) in the case of a particle on a 
worldline at constant radius. . 

we construct the mode 
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= 1 of the retarded Green function, 



Gret 
I 



{r,r';t), in the scalar case s — 



= r' = lOAf. We plot: (1) the BC contribution to G^'^*, i.e., Gf^, (2) the sum of Gf^ and the 



In Fig 
at the radii 

QNM contribution to G^*^* (taking into account the QNMs for the first 24 overtones), and (3) the full 'exact' G^'^*. 
We calculate Gf'~'' by integrating over the frequency the BC modes AGg obtained using two different methods in the 
'mid '-frequency and large-frequency regimes - in this case, the crossover frequency is at i> = 50 . In the 'mid'-frequency 
regime, we interpolate the values of the BC modes AGi obtained as described in the previous section. In the large- 
frequency regime we use the asymptotics of [5] for the BC modes. We calculate the QNM contribution to G^'^* using 
the method in [5S]. We obtain the 'exact' G^"^* by numerically integrati ng th e (l-l-l)-dimensional partial differential 
equation {—df + d^^ — V) 4>£{r, t) — (with the potential V given by Eq.(2.4)) for the £-mode 0£ of the field using the 
following initial data. We choose zero data for the inital value of the £-niode of the field. For the initial value of the 
time-derivative of the £-mode of the field, on the other hand, we choose a Gaussian distribution in 'peaked' at a 
certain value r^,o. From the Kirchhoff integral representation for the field (e.g., [T]), the solution (t)e{r, t) thus obtained 
should approximate the ^-mode of the retarded Green function, Gf'^*{r, tq; t), where vq = r(r*o) . We used a Gaussian 
width of approximately 0.2 A/ and we checked that the change in the numerical solution obtained by using smaller 
values of the width was negligible for our purposes. For the numerical integration of the (l+l)-dimensional partial 
differential equation we used Wardell's C-code available in |1D]. A slightly different version of this numerical approach 
using a Gaussian distribution (though using it as the source, rather than as initial data) has recently been successfully 
applied in [?T] in the full (3-|-l)-dimensional case. We observe from Fig 17 that the BC contribution becomes most 
significant for small values of the 'time' T = t — \r^, \ — \r'^, \ but, in the regime plotted, the BC contribution is always 
subdominant to the QNM contribution. The matching between the numerical solution and the sum of Gf^ plus QNM 
series is excellent. For T < neither the QNM series nor the BC integral is expected to converge separately [9]. 
Let us now define the '^-mode of the partial field' as 



/partial 



r) 



dt Gl^\r, 



r;t) 



(9.2) 



2\r,\ 



The contribution to the radial component of the self-force per unit charge in the case of a particle at constant radius 
from the £-mode of the Green function from the segment of the worldline lying between t = 2|r*| and t — oo is then 
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FIG. 16. Log-plot of '\2uq/{MW'^)\' (this is the branch cut mode AGe of Eq. (2.101 but without the fi factors) as a function 
of u. The vertical lines are located at the values of (minus) the imaginary part of the QNM frequencies - obtained from |3H] 
for the cases s = 1 and s — 2 and using [3^ for s — 0. The Wronskian W and the BC 'strength' q{i') have been calculated by 
evaluating the radial solutions at, respectively, r = 2.8A/ and r — 5M. (a) s — 0, £ = 1. (b) s = 1, ^ = 1. (c) s = 2, £ = 2: in 



this case we log-plot '|2i^g(!/ — fAs) M/W \\ i.e, we include an extra factor 
ai V = VAs)- 



VAS 



to account for the double zero of |W^| 
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obtained as: 



— = 2 — Pe{cos-f) — 



q 



dt 



^partial 

dr 



j^partial 



(9.3) 



This clearly only yields a partial contribution to the self-force from the ^-mode since we are integrating from t = 2|r»| 
instead of from t = 0+, as required in order to obtain the self- force. Because of the divergence of the BC and 
QNM contributions for T < 0, in order to obtain the contribution from the worldline segment for t : ^ 2|r*| 
we require a different method for calculating the Green function, such as a quasi-lo cal series (see, e.g., [42]). The 
BC contribution to i/)^'"^*'" is obtained by inserting Gf^ in the place of G^'^* in Eq.(9.2). In Fig 18 'a) we plot this 
contribution and its -derivative (evaluated using a central difference scheme) as functions of the radius in the case 
s = and £ = 1. In Fig 18 h) we plot the corresponding BC contribution to ff°,!'^^"'\r) / q in the static case 7 
dr/dt (1 - 2M/r)i/2. 



and 




FIG. 17. Mode £ = 1 of the retarded Green function, i.e., GY^{r,r' -,1), for s = and r — r' = lOM as a function of the 
'time' T = t — |r, | — Continuous blue curve: 'exact' numerical solution GY'^ obtained using [40]. Dashed black curve: 

). Dot-dashed red curve (overlapping with the continuous blue curve): Gf'~^ plus the corresponding QNM series 
contribution (summing overtone numbers n : — > 23) obtained using the method in [39] . The matching between the numerical 
solution and the sum of Gf ^ plus QNM series is excellent. 



Gf- of Eq.(2.9 




-0.00& 



FIG. 18. 'Partial' BC contribution to the ^-mode for the radial component of the self- force for s = Q a nd (. = 1 as a function 
of r«. (a) Continuous blue curve: BC contribut ion to the '^-mode of the partial field' (j)'^"'^^^"'' / M of Eq.(9.2l. Continuous red 
curve: its r, -derivative, (b) }i'^J^^'^\'r) / q of Eq.(9.3l in the static case (i.e., 7 = and dr/dt = (1 — 2M /r)^^). 
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X. DISCUSSION 



In this paper we have presented the first analytic method for calculating the branch cut modes in the non- asymptotic, 
'mid'-frequency regime in Schwarzschild space-time for fields of any integral spin. We have investigated their prop- 
erties, in particular regarding their relation to quasinormal mode frequencies and around the algebraically special 
frequency. We have applied our calculation of the BC modes to investigate their partial (i.e., from t > 2|r*|) contri- 
bution for one ^-mode to the self-force on a scalar charge moving on Schwarzschild background at constant radius. 
We have found that, for the particular case investigated, the BC contribution becomes larger as t approaches 2|r*| 
(where the high-frequency asymptotics of the BC modes become important) but the QNM contribution dominates 
the self-force at most times. 

In [43] the first successful application of the so-called method of matched expansions for the calculation of the 
self-force was achieved. This method consists on calculating the Green function at 'early times' using a quasilocal 
expansion and in the 'distant past' using a Fourier mode and multipole decomposition of the Green function as in 



Eq.(2.1 ). In [43] the method was applied to the specific case of a black hole toy model space-time, namely the Nariai 
space-time, where the Green function possesses QNMs but not a BC and so the Green function in the 'distant past' 
is fully determined by the QNM series. In Schwarzschild space-time, on the other hand, the QNM series must be 
complemented by a BC integral. In [li] the QNM series was calculated for large--^ and it was shown that it yields an 
interesting global singularity structure of the Green function. In this paper we have presented a method for calculating 
the BC integral and we have applied it to one £-mode. In [45] we plan to apply a calculation of the BC integral for 
all ^-modes, add it to a similar calculation of the QNM series and supplement it with a quasilocal expansion at 'early 
times' in order to calculate the full self-force in Schwarzschild space-time using the method of matched expansions. 

Another situation where it is important to investigate the contribution of the BC is that of the response of a black 
hole to an initial perturbation. In this case, we expect that the BC contributes at 'early' times (that is, for t close 
to |r*| + |r*|) as well as at late times (where a logarithmic behaviour precedes the known power-tail decay [6]). We 
investigate the latter in depth in [71. 

On the quantum side, the asymptotically constant spacing in the imaginary part of the highly-damped QNM 
frequencies led to suggestions of a link with the quantization of the black hole area [TOl |46]. In [9] we showed that, 
in the large-j> regime, the spacing in the imaginary part of the QNM frequencies asymptotically equals that of the 
zeros of the BC modes for all spins s = 0, 1 and 2. In this paper we have shown that, in the 'mid'-frequency regime, 
these two spacings also remain very close in the cases s = 1 and 2 studied, while they differ more significantly in the 
case s = 0. Intriguingly, the algebraically special frequency for gravitational perturbations plays a special role in the 
connection between QNMs and the BC, not only harbouring in its neighbourhood an almost purely imaginary QNM 
frequency but also marking the onset of the highly-damped regime for QNMs. In a different work, highly-damped 
QNMs in Kerr space-time have been interpreted as semiclassical bound states along a specific contour in the complex-r 
plane and have been linked to Hawking radiation [TT]. The least-damped QNMs have also been linked to Hawking 
radiation [17]. Given that QNM frequencies, particularly in the spin-1 case, 'approach' the branch cut 9] in the 
high-damping limit and that a connection with the branch cut also appears to exist in the 'mid'-frequency regime, it 
would be interesting to investigate whether branch cut modes may play any role in the quantum properties of black 
holes. 

An impending generalization of our current results in Schwarzschild is that to a rotating, Kerr black hole space- 
time. In principle, our method is readily generalizable to the rotating case, since Leaver's series representations for 
the radial solutions are already valid in Kerr [24] . Immediately, however, some significant differences appear with 



respect to the non-rotating Schwarzschild case. For example, the corresponding symmetry (2.8) in Kerr also involves 
a change in the sign of the azimuthal angular number, on which the radial solutions depend in the rotating case. As a 
consequence, the radial solution fi is not necessarily real along the branch cut and the discontinuity of across the 
branch cut is generally not only in its imaginary part but also in its real part, and so the BC 'strength' q(v) is not 
necessarily real-valued. To further spice up the analysis in Kerr, the angular eigenvalue has various branch points in 
the complex- frequency plane (see, e.g., [48]). This intricate and delicate structure of the modes in Kerr has various 
physical manifestations. We plan to investigate these issues in a future publication. 



Appendix A: Irregular confluent hypergeometric [/-function 



In this appendix we give some properties of the irregular confluent hypergeometric [/-function, which we have used 
in the main body of the paper. 
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A useful integral representation of the t/-function is given in, e.g., Eq. 13.4.4 [27] : 

1 r°° 

U(a,b,z)^—- dt e-''H''-\l + t)''-''-\ Re(a) > 0, |ph(z)| < 7r/2 (Al) 
r(a) Jo 

We require asymptotics for large values of the first argument of the [/-function. For large, negative values of the 
first argument we use Eqs.13.8.10 and 10.17.3 [27] to obtain 

—^1=] cos ^V2a?(6^^2^ + 7r ^a- ^ + , a -> -oo, 5 > 1, x > 

(A2) 

For large, positive values of the first argument we use Eqs. 13.8.8 and 10.25.3 j27j to obtain 

/ I \ l-b 



Uia,b,x)^ ^{xr-^e^^/'] -2Vk|ae'-^/ ^ a^+cx), 6<1, ^^arg(a;)e(-7r,+7r] (A3) 



Finally, from Eqs. 13.2.4 and 13.2.41 [57] we obtain the following expression for the discontinuity across the branch 
cut of the U -function, 

U{a, b, - z''U{a, 6, z) = ^.^f'"^^^^^" ,M {b - a, b, -ze^^^ (A4) 

^ ' 1 [a)L [1 + a — b) 
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